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Renormalization constants can be computed by means of Numerical Stochastic Perturbation Theory to 
two/three loops in lattice perturbation theory, both in the quenched approximation and in the full (unquenched) 
theory. As a case of study we report on the computation of renormalization constants of the propagator for Wilson 
fermions. We present our unquenched (Nf — 2) computations and compare the results with non perturbative 
determinations. 



INTRODUCTION & MOTIVATION 

Numerical Stochastic Perturbation Theory 
(NSPT) is a numerical method to perform 
Lattice Perturbation calculations. It has been 
shown to be very powerful in the quenched ap- 
proximation, but the technique is quite general 
and can be in principle applied to any theory 2 . 
In particular, it can be generalized to actually 
carry on unquenched QCD calculations 3 . The 
aim is to use NSPT to look at interesting quan- 
tities like improvement and renormalization co- 
efficients (quark propagator, currents, etc.), up 
to 2-3 loops and potentially for any fermionic ac- 
tion. In these pages we will present our results 
about the critical mass and the field renormaliza- 
tion constant (up to 3 and 2 loops respectively). 

1. COMPUTATIONAL SETUP 

For a comprehensive account of the strategy of 
unquenched NSPT we refer the reader to [3] and 
references therein. 

Simulations are performed on a APEmille crate 
(128 FPUs for 64 GFlops of peak performance), 
for Nf = 0, 2, 3, massless quarks Wilson action 
on a 32 4 lattice; gauge is fixed to the Landau 
condition using a Fourier accelerated algorithm. 
Stability with respect to single/double precision 



has been tested, as well as finite size effects are 
checked looking at smaller lattices on a PC clus- 
ter. So far we reached quite a wide configurations 
database (some hundreds, 1.8GB each) ready for 
any other observable. 

2. RESULTS 

Both critical mass and field renormalization 
constant are extracted from measurements of the 
quark propagator. Actually one averages S on 
the configurations and then gets T 2 by inversion. 
The structure of the result is as follow 



T 2 {p 2 ,m) 



where 



ip 1 — m — S(p 2 , to) 



(1) 
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£(p 2 , to) = S c + to E s (p 2 , to) + irf T, v (p 2 , m) (2) 

So, projecting our numerical results onto the 
7-matrices one recovers the last term in the sum 
and, from it, one can reach the field renormal- 
ization constant; in the same way, projeting onto 
the 7-identity one recovers the other two pieces 
of (0), which in (our) massless case reduce just to 
E c , the so-called critical mass: it is the additive 
renormalization for the mass one has to face for 
having broken chiral symmetry on the lattice. 
We present results only for the Nf — 2 case. 



1 



2 



-0.4 



-0.8 



-1.2 



-1.6 



2 9 

a p 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

Figure 1. Renormalized critical mass, (upper 
curve) and f}~ 2 coefficients; true data are repre- 
sented by the small points with errorbars; cir- 
cles represent the same points as they are recon- 
structed by the fitting procedure. 

2.1. Critical mass 

Perturbative expansion for E c is already ana- 
lytically known up to 2 loops | lloj . If we write 



e c = s«/r 



£< 3 >/r 



(3) 

,(2) _ 



for N f = 2 we have Y,c = 2.6057 and E^ 
4.293. In order to go beyond these coefficients 

(3) 

and get Ec , one must plug them in as countert- 
erms in the simulation because they contribute 
to higher loops results. So, if one looks at these 
terms in our results, one actually sees the addi- 
tive renormalized coefficients. This is just what 
is shown in Fig. ^ The non-smoothness of these 
curves is a consequence of having the continuous 
euclidean 0(4) symmetry broken on the lattice: 
at finite lattice spacing E is not simply a function 
of a 2 p 2 but it takes also corrections from other 
lattice invariants one can construct as powers of 
ap. Taking into account both the leading a 2 p 2 
ans higher order invariants, one can extrapolate 
the a — > limit value one is interested in. These 
values are just zero because of the counterterms 
subtraction. 

Fig. shows the same result, but for the 3-loop 
term. For it there are no counterterms plugged in 
the simulations (this is precisely the result one is 
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Figure 2. Critical mass, f3~ 3 coefficient; circles 
and points have the same meaning as in Fig. ^ 
the cross represents the extrapolated value. 

looking for), and so the extrapolation is not zero. 
The cross represents the extrapolated (a — + 0) 
value, which reads 



E^ = 11.79 



(2) 
(5) 



(4) 



in which the errors are estimated by looking at 
the stability of the fitted result with respect to 
varying the number of points included in the fit 
and the higher orders lattice invariants taken into 
account. 

2.2. Field renormalization constant 

As said, the field renormalization constant 
Z q can be extracted looking at 5*" 1 along the 
7-matrices 

_ i Tr^S- 1 ) 
q ~ 4 p 2 

= 1 + Z^p- 1 + Z^p- 2 + . . . (5) 

In contrast with the case of the critical mass, in 
this case we have to face logarithmic divergences, 
log(a 2 p 2 ), coming from Renormalization Group 
Theory. These divergences are ruled by anoma- 
lous dimensions and so, because we are in Landau 
gauge, we have just no log-terms at 1-loop and 
only a simple log-term at 2-loops (i.e. without 
log 2 -term). 

So, at 1-loop one has no divergences to subtract 
and one can extrapolate directly to the a — > 
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Figure 3. Field renormalization constant, /3 1 
coefficient; marks have the same meanings as in 
Fig. El 

limit. The result we found is shown in Fig. [3] and 

reads 



ZjU = -0.8432]. 



(6) 



This result is in perfect agreement with the known 
value Z ( q ] = -0.843... 

At 2-loops one has to subtract the known 
log-term before taking the continuum limit. In 
Fig. 0] we plot both the crude data and the data 
with the log divergence subtracted: the extrapo- 
lation for the latter reads 

Zf> = -1.34(2). (7) 

We can compare this new result with a non- 
perturbative value 6 , which at /3 = 5.8 and 
a/i = 1 reads ~ P) = 0.786(5). Our 2-loop ex- 
pansion add up to Z\ ' — 0.815. So the discrep- 
ancy reduces from 9% (1-loop) to 3-4%. In order 
to improve convergence, one can switch from a 
bare to a better-converging coupling expansion 
(such as in the so-called Boosted Perturbation 
Theory). On the other hand we already have the 
signal for the 3- loop coefficient. 

3. CONCLUSION & PERSPECTIVES 

We are currently refining the statistic in order 
to reduce statistical errors and to see the 3-loop 
coefficient for Z„. 
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Figure 4. Field renormalization constant, /3~ 2 
coefficient; marks have the same meanings as in 
Fig. |2 crude data are on the lower curve; on the 
upper curve the log-divergence is subtracted. 

Just the same strategy used here is needed to 
measure fermion currents and the renormaliza- 
tion constants related to them (and we are actu- 
ally working at it) as well as to investigate the 
twist-mass QCD action. 

Much the same strategy can also give improve- 
ment coefficients. 

Almost the same strategy may also allow to 
perform (quenched) overlap fermion calculations 
(and we are actually working on it as well). 
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